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ABSTRACT 

We present the results of a study for galaxy orbits in galaxy clusters using 
a spectroscopic sample of galaxies in Sloan Digital Sky Survey (SDSS) and 2dF 
Galaxy Redshift Survey (2dFGRS). We have determined the member galaxies of 
Abell clusters covered by these surveys using the galaxies' redshift and positional 
data. We have selected 10 clusters using three criteria: the number of member 
galaxies is greater than or equal to 40, the spatial coverage is complete, and X- 
ray mass profile is available in the literature. We derive the radial profile of the 
galaxy number density and velocity dispersion using all, early-type, and late-type 
galaxies for each cluster. We have investigated the galaxy orbits for our sample 
clusters with constant and variable velocity anisotropics over the clustercentric 
distance using Jeans equation. Using all member galaxies, the galaxy orbits are 
found to be isotropic within the uncertainty for most of sample clusters, although 
it is difficult to conclude strongly for some clusters due the large errors and the 
variation as a function of the clustercentric distance in the calculated velocity 
anisotropics. We investigated the orbital difference between early-type and late- 
type galaxies for four sample clusters, and found no significant difference between 
them. 

Subject headings: galaxies: clusters: general - galaxies: clusters: individual 
(A85, A779, A1650, A1651, A1795, A1800, A2034, A2199, A2670, and A2734) - 
galaxies: kinematics and dynamics 
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Introduction 



The mass estimate of a galaxy cluste r, the la r gest K ravitationally bounded system in the 
universe, was given for the first time by IZwickyI (jl933l ). which showed an indication of the 
existence of dark matter. In general, the cluster mass ranges from IO^^Mq for small groups 
to lO^^MfT) for rich clusters, which is in the form of gala xies (1 — 2% of the to tal mass, e.g., 
Lin et al.ll2003l ). hot X-ray emitting gas (5 — 15%, e.g., IVikhlinin et al.ll2006l ). dark matter 
(80 — 90%) that affects the galaxy cluster through the gravitation only. 

To determine the amount and distribution of underlying dark matter in clusters, Jeans 
analysis using the positional and velocity data of cluster galaxies has been usually adopted 
among several mass estimation methods (iKatgert et al.ll2004l : lBiviandl2006l ). However, pre- 
vious studies had to assume the galaxy orbit in prior to derive the mass profile of galaxy 
clusters , due to the 'mass-orbit' degeneracy in the velocity dispersion profile (VDP) (e.g., 
MerrittI 119871 ). This degeneracy can be broken by an orbit analysis with an independent 
mass determination based on X-ray or lensing studies, or by an anal ysis of higher order mo- 



ments of velocity distribution : velocity dispersion and kurtosis (e.g..lLokas fc Mamon 



Lokas et al.ll2006l ). Gauss- Hermite moments (e.g.. Ivan der Marel et al. 
2004[ ). 



2000 



2003 



Katgert et al. 



The analysis of ga l axy or bits for clusters was made for the Coma cluster for the first 
time by iKent fc GunnI (119821 ) who used about 300 galaxy velocities. They showed that 
the galaxy orbits in Coma can not be primarily radial, and a significant fraction of kinetic 
energy must be in tangential direction even at large radii. With an aid of large redshift sur- 
veys for cluster galaxies such as Canadian Network for Obs ervational Cosmolo gy (CNOC; 
Yee et al.lE996[ ). ESQ Nearby Abell Cluster Surve y (ENACS; iKateert et al.lll996h . and Clus- 
ter and Infall Region Nearby Survey f CAIR NS: iSines et aljE)03h. extensive analysis for 
galaxy orbits in clusters has been performed. ICarlberg et al. (1l997al lb[) analyzed 14 CNOC 
cluster s at z = 0.17 — 0. 5 4, an d found that the galaxy orbits are isotropic or modest ra- 



dial. Ivan der Marel et al.l (120001). using t h e sara e CNOC data, concluded that their best-fit 



model is close to isotropic. iRines et al.l (120031 ) secured 15,000 galaxies in CAIRNS data, 
and concluded that the galaxy orbits in their ensemble cluster are also consistent with being 
isotropic. 

On the other hand, the difference in galaxy orbits and in galaxy velocity dispersions 
am ong galaxy types provide important c lu es to understand the formation history of the clus- 
ter (IBiviano fc Katgertll2004j : lGotoll2005l ). iBiviano et al.l (119971 ) studied kinematic difference 
between emission-line galaxies (ELGs) and non-ELGs using ENACS data, and f ound that the 



VDP of ELGs is consistent with being radial, while that of non-ELGs is not. lAdami et al. 



(1l998l ). based on simple modelling of observed VDP for about 2000 galaxies in 40 regular 
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clusters, reported that the orbits of elhpticals are mostly tangential in the c luster core and 



are ne arly isotropic outside, while those of spirals are predominantly radial. iMahdavi et al. 



( 119991 ) analyzed a sample of 20 galaxy groups from Center for Astrophysics redshift survey. 
They found that star-forming galaxies or ELGs have moderately radi al orbits, while old 



or abs orption-line galaxies isotropic orbits within the errors. Recently, iBiviano fc Katgert 



( I2OO4J ). using the EN ACS data, studied the galaxy orbits of all galaxy classes (the brightest 
ellipticals, other ellipticals together with SO, early-type spirals, late-type spirals, and irregu- 
lars) for the first time. They reported that the brightest ellipticals do not yie ld equilibrium 



solution, other ellipticals together with SO have isotropic orbits as given in iKatgert et al. 



( 120041 ). and early spirals are consiste nt with isotropic orbits, whi le late spirals prefer radial 
orbits to isotropic orbits. In contrast. iRamirez &: de Souzal (Il998l ). using nearby {z < 0.055) 
Abell clusters, concluded that the orbits of elliptical galaxies in clusters are close to radial, 
while those of spirals more tangent ial or isotropic. For nine interme diate jz ~ 0.1 — 0.4) 



CNOC cluster. iRamirez et al.l (I2OOOI ) obtained similar results to those of lRamirez fc de Souza 
(119981 ): bulge-dominated galaxies have more eccentric orbits than disk-dominated galaxies 
do. The cause for t he different results o f Ramirez et al. co r npare d with other studies is 
discussed in detail in Ivan der Marel et al.l (I2OOOI ) and iBivianol (120021 ). 



Above results of galaxy orbits are based on a composite clusters prepared by combining 
data for disparate clusters that might have different formation histories. Although it is 
helpful to make the composite cluster in order to overcome problems of 1) a limited number 
of measured galaxy velocities per cluster and 2) an application of spherical Jeans equation 
to an asymmetric cluster, the composite cluster might be significantly different from any 
real cluster. In addition, previous studies were not based on the independent measurements 
of cluster mass profiles, but based o n the optical galaxy data. These problems can be 



overcome. 



Natarajan fc KneibI (119961 ) demonstrated the galaxy orbits can be constrained 
using an independently determined mass profile. For A2218, they used mass estimates 
derived by c ombining stroiig and weak lensing effects, and found evidence for an anisotropic 
core. Later, iBenatov et al.l (120061 ) extended the study of galaxy orbits to five clusters using 
mass estimates from X-ray data. They reported that the galaxies in five clusters have diverse 
orbits, and the orbital profiles in observed and simulated clusters appear to be different: 
simulated clusters show preferentially tangential orbits. However, they did not divide the 
galaxy sample into different types to investigate the orbital properties of different galaxy 
types. 



Recently, galaxy redshift surveys such as the Sloan Di gital Sky Suryey ('S DSS: lYork et al. 



2000) and the 2dF Galaxy Redshift Survey (2dFGRS; ICoUess et al.ll200lh have provided 
redshift data for large sa mples of galaxies. In addition, X-ray mass profiles for numerous 



clusters become available (iReiprich fc Bohringerll2002l : ISanderson et al.ll2003l : iDemarco et al. 
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20031 : iBrownstein fc MoffatlboOQ : IVikhlinin et al.lbood : IVoigt fc Fabianll2006r ) . By identifying 



member galaxies in clusters using these redshift and positional data and adopting X-ray mass 
profiles, we can investigate the galaxy orbits in clusters for a large sample of galaxy clusters 
without making the composite cluster. We can also investigate the orbital properties for 
different galaxy types. 

In this paper, we present the results of a study for the orbits of galaxies in galaxy clusters, 
using a spectroscopic sample of galaxies in SDSS and 2dFGRS and X-ray mass profiles in 
the literature. Section [2] describes the sample of galaxies and clusters used in this study. 
Radial profiles of galaxy number density and velocity dispersion are derived in ^ Analysis 
of the galaxy orbits and global kinematics are given in §1] and ^ respectively. Discussion 
and summary are given in ^ and §3 respectively. Throughout, we adopt cosmological 
parameters h = 0.7, Q\ = 0.7, and Qm = 0.3. 



Data 



2.1. Galaxy Sample 



We used data fr om the spectroscop ic sample of galaxies in the Legacy surv ey of SDSS 
Sixth Data Releas^ ( Strauss et al.|[2007 ). and in the 2dFGRS Final Data Release ( CoUess et al 
20011 ). SDSS is one of the largest imag ing and spectroscopic survey, which cover s more than 



a qua rter of the sky (lYork et al.ll2000l ) using a dedicated, 2.5-meter telescope ( iGunn et al. 



20061 ) at Apache Point Observatory, New Mexico. The telescope is equipped with a wide- 
field dri ft-scanning mosaic CCD camera that can image 1.5 deg^ of sky at a time for imaging 
survey (iGunn et al.lll998l ). and a pair of fiber-fed spectrographs that can measure spectra 



more than 600 objects in a single observ ation with 3° field of view (jUomoto et al.l Il999 



Castander et al 



2001 



Blanton et al.ll2003r). Extensive d escription of SDSS data products is 



given by lYork et al.l (120001 ) and IStoughton et al.l (120021 ) , and we only give a brief summary 



of the SDSS data. The imaging survey is carried out with fiv e broad bands {ugriz) centered 
at 3551, 4686, 6166, 7480, and 8932 A jFukurita et al.lll996h . The magnitude limits of five 
bands for point sources are 22.0, 22.2, 22.2, 21.3, and 20.5, respectively. The calibration 
of data p roducts from imagi ng survey is given by sever al authors: data r eduction and pho- 



Ivezic et al. 



2004 



tometry (ILupton et al.l 120021). photom etric calibration (jHogg et al.ll200ll: ISmith et al.l 12002 



Tucker et al.l l2006l ) . and astrometric calibration (jPier et al. 



200 



The 

spectroscopic survey obtains spectra covering 3800—9200 A, with a wavelength resolution. 



^Access to Data Release 6 can be found on the SDSS Web site ( jhttp: / /www. sdss.org/dr6 1 . 
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A/5A ~ 1800. There are three kinds of spectroscopic targets: main galaxy sample, lumi- 
nous red galaxies (LRG), and quasar. The main galaxy sample consists of the galaxies with 
r-band Petro sian magnitude rppt, < 17.77, which are corrected for Galactic foreground red- 
dening using ISchlegel et al.l (119981 ) . The galaxies with r-band Petrosian half-light surface 
brightness /X50 > 24.5 mag arcsec"^ are rejected; the number ratio of rejected galaxies to all 
galaxies is only ~ 1%. In addition, the galaxies with 3" fiber magnitudes brighter than 15 
magnitude in g or r, or brighter than 14.5 magnitude in i are rejected to avoid the problems 
of saturation and cros s-talk in the spectro graphs. A detailed description for main galaxy 
sample is presented in IStrauss et al. fl2002h . 



We used the main galaxy sample of 528,000 galaxies with measured velocities for this 
study. The median redshift for this spectroscopic sample is 0.11. The uncertainties of the 
redshift measurements are ~30 km s^^. A redshift confidence parameter (zConf) is assigned 
from to 1 in the SDSS catalogs. We used only those galaxies with zConf > 0.65. 

2dFGRS is a spectroscopic survey for near ly 246,000 ga.laxies selected in the photo- 
graphic 6j band from the APM galaxy catalog (jColless et al.ll200ll ). The survey uses the 
Two-degree Field (2dF) multifibre spectrograph on the Anglo- Australian Teles cope, which 



meas ure spectra about 400 objects in a single observation with 2° field of view (iLewis et al. 



20021 ). The wavelength coverage of the spectra is 3600—8000 A with a resolution of 9 A. The 
spectroscopic targets are the galaxies with extinction corrected magnitude 6j < 19.45. The 
survey coverage consists of two delination stripes and 100 random fields, which cover in total 
over 2000 deg^. The median redshift for the sample of galaxies is similar to that of SDSS. A 
redshift quality parameter Q is assigned in the range 1—5. The redshift measurements with 
Q>3 are 98.4% reliable, and have an overall rms uncertainty of ~85 km s~^. We used only 
the galaxies with Q>3 in this study. 

Since some galaxies are covered in both SDSS and 2dFGRS, we matched the galaxies 
found in SDSS with those in 2dFGRS to make a master catalog. The mean difference in 
radial velocity At> (= fsDss — ^^2dFGRs) between the SDSS and 2dFGRS measurements for the 
29,200 matched galaxies is estimated to be ~ 13 km s~^. We corrected the SDSS velocities 
by this mean difference, and used the average value of the velocities measured in SDSS and 
in 2dFGRS for further analysis. 

To investigate the difference of galaxy orbits among subsamples, we divide the galax- 
ies into early-type and late-type galaxies using spectroscopic parameters provided by each 
survey. For 2dFGRS galaxies, r] parameter, which is a linear combination of emission and 
absorption components within t he spectrum derived fr om the Principal Component Anal- 
ysis, is assigned to each galaxy (IMadgwick et al.ll2002l ). This para meter denotes the ratio 
of the present to the past star-formation activity in each galaxy. iMadgwick et al.l (|2002| ) 
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showe d that the r? par ameter correlates well with galaxy morphology using the Kennicutt 
Atlas (iKennicuttl Il992l ) (see their Fig. 4). They reported that the galaxies with rj < —1.4 
are corresponding to E/SO and Sa and those with 77 > — 1.4 to Sb and Scd. Similarly, for 
SDSS galaxies, eclass parameter, which is a projection of the first three principal compo- 
nents of the spectrum, is assigned to each galaxy. This parameter has values from about 
—0.35 (correspondi ng to early-type galaxies) t o 0.5 (late-types), which was used to classify 
SDSS galaxies (e.g.. iBernardi et al.ll2003l . l2006l ). In Figured! we plot the eclass parameter 
in SDSS versus the r] parameter in 2dFGRS in order to determine the eclass value corre- 
sponding to the division value of r/ = — 1.4. It is seen that two parameters correlate well and 
the majority of galaxies are located around [rj, eclass) = (— 2.65,— 0.135). We determined the 
relation equation between two parameters using the ordinary least-square bisector method 
(llsobe et al.lll990l ). The fit was done for the galaxies with rj < 1.1 and eclass< 0.09. 



eclasssDss = 0.0631(±0.0004) r/2dFGRS + 0.0244(±0.0008) 



(1) 



This relation indicates that the division value of 77 = —1.4 in 2dFGRS is equivalent 
to that of eclass=— 0.0640 in SDSS. Therefore, we classified the galaxies into early-type 
galaxies if (i) rj < —1.4 when the galaxy is surveyed only in 2dFGRS, (ii) eclass< —0.0640 
when the galaxy is surveyed only in SDSS, and (iii) eclass< —0.0640 and 77 < —1.4 when 
the galaxy is surveyed both in SDSS and 2dFGRS. The rest were classified into late-type 
galaxies. Some galaxies without estimation of rj parameter (77 = —99.9) were not used for 
the analysis of subsamples. 



2.2. Cluster Sample and Galaxy Membership in Clusters 



We used the Abell catalog of galaxy clusters (lAbell et al.lll989l ) to identify cluster galax- 
ies in the survey data. Among the Abell clusters, we selected those that have known spec- 
troscopic redshifts in the NASA/IPAC Extragalactic Database (NED). Finally, we selected 
731 and 230 clusters located within the survey regions of SDSS and 2dFGRS, respectively, 
as a sample of clusters for further analysis. 

In order to de termine the r nemb ership of galaxies in a cluster, we used the "shifting 
gapper" method of iFadda et al.l (119961 ) as used also for the study of global rotation of galaxy 
clusters ( Hwang fc Leel 2007 ). In the plot of radial velocity versus clustercentric distance of 
galaxies for a given cluster, we selected the member galaxies using a velocity gap of 950 km 
s~^ and a distance bin of 0.2 Mpc shifting along the distance from the cluster center. We used 
a larger bin width if the number of galaxies in a bin was less than 15. We applied this method 
to the galaxies within the radius at which the distance between adjacent galaxies is larger 
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than 0.1 Mpc. If there are no adjacent galaxies with > 0.1 Mpc, we stopped the procedure 
at the radius of 3.5 Mpc. We iterated the procedure until the number of cluster members is 
converged. Finally, we selected 113 galaxy clusters in which the number of member galaxies 
is greater than or equal to 40 for further analysis. 

Since it is necessary to determine the galaxy orbits using cluster mass profiles from X-ray 
data, we selected 21 galaxy clusters whose X-ray mass profiles are available in the literature 



Reiprich fc Bohringe 1I I2OO2I : ISanderson et ahlbood : IPemarco et al. I I2OO3I : iBrownstein fc MoffatI 



20061). Then, we rejected 9 clusters (A1775, 2052, 2063, 2142, 2147, 2244, 2255, 4038, and 
4059) that appear to be in the stage of interacting or merging in the plot of galaxy velocity 
versus clustercentric distance, and 2 clusters (A119 and A1656) whose survey coverages are 
not complete. Although A85 was partially surveyed, we included it into our sample since the 
uncovered region is only the small outer region of the entire cluster. Finally, we obtained a 
sample of 10 clusters that will be used for investigating the galaxy orbits. 

Table [1] lists our sample clusters with Abell identification, right ascension and decli- 
nation, Bautz- Morgan (B-M) type, the survey in which the clu ster is covered, the redshift 
derived in this study (the biweight location of iBeers et al.lll990l ). the phy sical extent corre- 
sponding to one arcmin, and the velocity dispersion (the biweight scale of iBeers et al.lll990l ) 
and the number of galaxies for all, early-type, and late-type galaxies. Our sample clusters 
are found from z = 0.023 to 2; = 0.113 and have 78—754 member galaxies. The number of 
early-types is usually larger than that of late-types in clusters, and eight clusters have B-M 
type of I or II. Figure [2] shows plots of radial velocity as a function of clustercentric distance 
of galaxies and the velocity distributions for the 10 clusters. In Figure [3], we show the spatial 
distribution of cluster galaxies with measured velocities. For most clusters, early-types are 
centrally concentrated, while late- types are not. 

Since the presence of substructure in clusters can affect the determination of the VDPs 
and the galaxy number density profiles, it is necessary to secure the galaxy sample outside 
the substructure in order to determine the galaxy orbits properly. Therefore, we show, for 
our sample cluster, the distribution of 6 that indicates local deviations from the systemic 
velocity (fsys) and dispersion ((Xp) of the entire cluster in Figure|l](see §5.2l for the explanation 
of 6 in detail). We selected the galaxies with 6 < 2.0 that are regarded as the galaxies in the 
cluster main body for the orbit analysis, and compared the results for the different choice of 
S value in ^6.2[ 
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3. Observed Properties of Galaxy Clusters 
3.1. X-ray Mass Profile 



We present X-ray mass profiles for our sample clusters in Figure O It appears that 
mass profiles from various references agree well as a whole in each cluster. In Table [2l 
we list the X-ray luminosity and the parameters for the X-ray mass profile. The first and 
second columns give Abell identification and the X-ray luminosity, respectively. The third 
and fourth columns are f3 parameter and core r adius (rr) derived from the standa rd /9-model 
of cluster gas density profile, respectively fe.g.. ICavaliere fc Fusco-Femiandll976r). The fifth 
and sixth columns are isothermal gas temperature ( Tx) i n iReiprich fc Bohringerl (120021 ) , and 
central gas temperature [T(0)] in lSanderson et al.l (120031 ). respectively. The seventh column 
gives the polytropic index (7) used in equation (jlj), and final column gives the reference 
of these parameters adopted in this study. Then, the gravitational cluster mass Mtot{r) is 
determined under the assumption of hydrostatic equilibrium through the equation. 



Mtot(r) 



kT{r)r f d Inp(r) d InT(r) 



Gfinip 



d Inr 



d Inr 



(2) 



where k is the Boltzmann constant, G is the gravitational constant, p (= 0.61) is the 
mean molecular weight, rrip is the proton mass. The gas density profile is given by. 



p(r) = p(0) 



(3) 



For the clusters in IReiprich fc Bohringerl ( 20od ). we use isothermal gas temperature (Tx) 
as T(r). For the clusters in lSanderson et al.l (l2003l ). we use non-isothermal gas temperature 
profile that is linked to the gas density via a polytropic equation of state. 



T(r) = T(0) 



1+ - 



/3(7-l) 



(4) 



where 7 is the polytropic in dex and and B are a s defined previously. For A2034, the X-ray 
mass profile is available in iDemarco et al.l (120031 ). However, they present a different form 
of gas density profile from the standard /5-model, so those parameters were not included in 
Table H 

It is worth noting that the X-ray mass profiles may not be well determined outside the 
outer significance radius of the cluster, ry, because the X-r ay source count drops below the 
Poissonian la error at r > (IReiprich fc Bohringerl l2002l ) . It is also noted that the Jeans 
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equation may not be applicable beyond a radius of r2oo (usually called virial radius), because 
the dynamical equilibrium of the cluster is not guaranteed at this region. Therefore, careful 
interpretation is needed beyond rx and r2oo- We draw the vertical lines at the radius of rx 
and r2oo in Figures 0, [TTl [T21 [13 [HI and [T3 The radius r2oo that contains an over density 
200pc wher e Pc is the critic a l densi ty of the Universe, is computed for each cluster using the 
equation in ICarlberg et al.l (Il997al ): 



200 



(5) 



where ad is a velocity dispersion of the cluster and the Hubble parameter at z is H^{z) = 
i^Q [f2jv/(H- -z)'^ + ^^fc(l + 2;)^ + ^^a] (|Peebleslll993l ). Qm, ^k, and are dimensionless density 
parameters. 



3.2. GalcLxy Number Density Profile 



We derived the galaxy number density profile for each cluster using the member galaxies 
with S < 2.0 selected in §2.21 Since SDSS and 2dFGRS are magnitude limit surveys (rpet < 
17.77 for SDSS and bj < 19.45 for 2dFGRS) and the spectroscopic sample of galaxies are 
nearly complete within the magnitude limit, we can derive the galaxy number density profile 
for each cluster using spectroscopically selected member galaxies. 

We display the galaxy number density profiles of all, early-type, and late-type galaxies 
for each cluste r in Figures M and [71 We de-projected the observed number density using 
the method iii [McLaughlin T1999L and fitt e d the de-projected density profile with those of 
Navarro et al.[ (Il997[ . NFW) and [Hernquist[ ([19901), which are represented by. 



(6) 



where t'g(r) is a three dimensional density profile of the cluster galaxies, is a scale 
radius, and a = 2 for the NFW profile an d a = 3 for t he He rnquist profile. We also fitted the 



de-projected density profile with that of [Kim et al.[ (l2005l KEK05), which is derived from 
the isotropic and isothermal galaxy orbits in clusters under the NFW distribution of dark 
matter. It is represented by. 



z/g(r) = 1^0 

where the dimensionless parameter rjg is defined by 

2^71^0 



(7) 



Vg 
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We used as a characteristic mass Mq, X-ray mass estimate shown in Figure [5] at a radius 
of r2oo, and used as a one-dimensional velocity dispersion a,., observed velocity dispersion cTp 
shown in Table [TJ 

The solid, dashed, and dot-dashed lines represent the projected best fit curves of the 
NFW, Hernquist, and KEK05 profiles, respectively. We did not fit the profiles that give 
unstable de-projection, and did not include them in Figures [6] and [71 For the clusters (A1650, 
1651, and 2670) surveyed in both SDSS and 2dFGRS, we used one data set (2dFGRS for 
A1650 and A1651, SDSS for A2670) for deriving number density profile since the spatial 
coverage of SDSS or 2dFGRS is not complete depending on clusters. For A85, the fit was 
done using only the galaxies in the inner region [R < 26') because of incomplete coverage 
(see Fig. ED. 

The fitting results for our sample are summarized in Table [31 It is found that the scale 
radius of late-types in all profiles is larger than that of the early-types, showing that the 
latter is more concentrated toward the cluster center than the former. The steeper density 
profile of earl y-types compar ed with late-types is a clear indication of the morphology- density 



relation (e.g.. lDresslerlll980l ). 



To test the effect of the incompleteness of spectroscopic sample on the galaxy number 
density profiles, we investigate the color magnitude diagram of the galaxies in the photomet- 
ric and spectroscopic samples. It is well known that optical colors of early-type galaxies are 
strongly correlate with absolute magnitudes (color-magnitude relation, CMR), in the sense 
that t he brighter galaxies are likely t o be redder than the fainter galaxies (lBaumlll959l : iFaber 
19731 : IVisvanathan fc Sandagelll977l ). In Figure [H we plot the g — r versus r for the SDSS 
cluster galaxies and the Bj — Rp versus Rp for the 2dFGRS cluster galaxies. The cluster 
galaxies within the Rmaxf^ from the cluster center selected in §2.21 are shown in company 
with the photometric sample of galaxies without measured velocities in the same region. 
Rmax is the largest clustercentric distance of the member galaxies, and varies depending on 
clusters. It shows that the colors correlate well with the observed magnitudes. The linear fits 
derived from repeated one sigma clipping using early-type galaxies selected in §2.11 are over- 
laid. The average of the CMR slopes for our sample clusters is found to be —0.019 ± 0.006 
for SDSS clusters and —0.037 ± 0.024 for 2dFGRS clusters, which is consistent with the 



results of preyious studies based on the same data (IDe Propris et al.l [20041 : [Hogg et al.ll2004 



Gallazzi et al.|[2006[ : [Aguerri et al.[[2007[ ). It appears that most bright (r < 14 mag) galaxies 
following the CMR in the photometric sample are selected as member galaxies, implying 
that the incompleteness of spectroscopic sample in the cluster center does not affect the 
determination of the projected galaxy number density in Figures M and [3 



3.3. Velocity Dispersion Profile 



We present, in Figures [9] and [TOl tlie VDPs for all, early-type, and late- type galaxies 
with 6 < 2.0 in each cluster. We computed the velocity dispersion of the galaxies lying 
within a bin with fixed radial width, AR = Rmax/^^ for all galaxies and AR = Rmax/^ 
for the early-type and late-type galaxies as increasing the bin center by a fixed step width, 
6R = Rmax/4:0. We stopped the calculation when the number of galaxies in a bin is less than 
9, and used larger bin size if the number of computed dispersions in a cluster is less than 
five. We do not present the profiles for the subsample of the cluster whose final number of 
computed dispersions is less than five. The velocity dispersions for our sample cluster are 
not similar, but show diverse features. In addition, they are not always constant throughout 
the radius. 



4. Galaxy Orbits in Clusters 

If we assume spherical symmetry of a collisionless galaxy cluster, we can apply the Jeans 
equation in the absence of rotation to the dynamical analysis of the cluster. The spherical 
Jeans equation is 



d , 



r)a,{r)) + 



tot\ 



(9) 



where arir) is the radial component of velocity dispersion, [3orhij') = 1 — (^lij) / (^ri''^) 
the velocity anisotropj]^, G is the gravitational con stant, and Mtntiv) is the to tal gravitating 
mass contained within a sphere of radius r (e.g. iBinney fc Tremaind 119871 ). (Te{f^) is the 
tangential component of velocity dispersion that is equal to the azimuthal component, cT0(r), 
i n the absence of clu ster rotation. However, cluster rotation is not negligible in some clusters 
jHwang fc LeelboOTh . 



With an aid of independent determination of the cluster mass profile using X-ray data, 
we determine the velocity anisotropy (orbits) of cluster galaxies using two methods: (1) from 
the c omparison of the calculated VDP us ing the Jeans equation with the measured VDP 



[e.g., 



Cote et al.l l200ll: iHwang et al 



Jeans equation (e.g.. iNatarajan fc KneibI 119961 : ISenatov et allboOfih . 



20071). and (2) to calculate d irectly /9orb('") using the 



^It is noted that /3orb('") is velocity anisotropy, while /3 is a parameter used for X-ray mass profile in Wi.W 
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4.1. Method 1 

Our strategy to determine the velocity anisotropy of cluster galaxies is as follows: (1) 
With the galaxy number density profile, z^g(r), of all, early-type, and late-type galaxies and 
the mass profile, Mtot{r), in hand, assuming constant velocity anisotropy, Porh, over the radius 
in prior, we compute the theoretical projected VDP, ap{R) and theoretical projected aperture 
VDP, crap(< R), using the Jeans equation; (2) From the comparison of these calculated VDPs 
with measured VDPs, we determine the velocity anisotropy of cluster galaxies. 

We begin by deriving the theoretical projected dispersion profiles. The equation ([9]), 
spherical Jeans equation can be solved for the radial component of velocity dispersion, ar{r): 



exp 



2A 



'orb 



dr 



Z/p. 



tot 



exp 



2A 



orb 



X 



dx I dx 



(10) 



Then the projected VDP, ap{R) can be derived by 



9 /"OO / 



r dr 



- i?2 



where R is the projected clustercentric distance, and Sg(-R) is the projected density 
profile that is a projection of the three-dimensional density profile t'g(r): 



Eg(i?) 



r dr 



R 



Vr2 - i?2 



(12) 



The projected aperture dispersion profile, crap(< R), which is the velocity dispersion of 
all objects interior to a given projected radius R, can be computed by 



R 



j:g{R')a'{R') R' dR' 



R„ 



R 



-^min 



Sg(i?') R' dR' 



-1 



(13) 



where -Rmin is the projected clustercentric distance of the innermost data point for the 
cluster galaxies. 

For all galaxies in a cluster, we present the measured VDP compared with the calculated 
VDP for different velocity anisotropics in Figures [11], [12], and [131 The top panels show the 
projected VDPs, and the middle panels show the projected aperture VDPs. The measured 
dispersion data taken from Figures [9] and [10] are shown by filled circles along with their 



- 13 - 



confidence intervals. The projected aperture dispersion profiles are plotted using a similar 
fashion to the case of the top panel. Although it is difficult to distinguish the velocity 
anisotropy clearly in the top panel (the middle panels show more stable result), it appears 
that the galaxy orbits in clusters show diverse patterns. Three results based on NFW, 
Hernquist, and KEK05 profiles of galaxy number density appear to be similar for a cluster. 
In Figures [H] and [15], we show the results of similar analysis for early-type and late-type 
galaxies, respectively. The difference between subsamples in a cluster is not clearly seen at 
this stage, and we discuss in detail the individual cluster in Appendix A. 



4.2. Method 2 

We combine the equation Q), spherical Jeans equation with the equation (ITTl) that 
defines the projected VDP in order to compute PoTh{i^) as a function of the radius directly. 

Replacing Porh in equation (ITTi) with that in equation we get. 



PC 

JR 



rvgaldr ^ f°°d{uga^) dr 



: dr 



R v/H^i?^ 2 



R 



dr - i?2 



(14) 



T his equation can be redu ced by integrating the second term on the right-hand side by 
parts (IBinney fc Mamonlll982l ). 



T.g{R)al{R) - I 

JR 



y 2 



R - R^ dr 



■ dr. 



(15) 



Then, Vgol can be expressed as a sum of four integrals (see iBicknell et al.l Il989l for 



detail). 



Vg{r)al = \{x)-\^{r)^\{r)-\{r), 



(16) 
(17) 
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2 



r 



Hr) = — I GMtot(r)z/,rdr, (1^ 







Hr) = ^ r Rj:,{R)al{R)dR, (19) 
r Jo 

2 f°° / r r \ 

Ur) = —j^ R^^i^Hi^i^/^T^ - (20) 

After computing z/gcx^, we finally obtain Porh{r) from the Jeans equation, 



^"^^^""^ " + — ('^f^r) • (21) 



Practically, all integrations up to infinity were performed using a large radius, Rt = 
1.5Rmax at which both i/gcr^ and Sg(i?) approach to zero. In addition, input VDPs derived 
in §3.31 were smoothed for stable computation. We compute Porh using the above equation, 
and present the results in the bottom panels of Figures [T H [T2 | and [131 The errors of /5orb are 
computed using the upper and lower confidence intervals of the VDPs shown in Figures [9] and 
[TUl The behavior of /5orb computed in this Section is similar to that in §4.11 It is noted that 
the calculated /3orb may not be reliable in the very inner region where the observed VDPs 
are not available (interpolation of the observed VDPs in the outer region is used) and X-ray 
mass profiles are not properly determined due to the low resolution of X-ray instrument. 
Three results based on NFW, Hernquist, and KEK05 profiles of galaxy number density are 
also similar for a cluster. We also show the results of similar analysis for early-type and 
late-type galaxies in Figures [H] and [151 respectively. Detailed discussion of the individual 
cluster is given in Appendix A. 



5. Global Cluster Properties 

5.1. Cluster Morphology 

In order to investigate the connection between cluster morphology and cluster dynamics 
it is useful to determine the ellipticity and orientation of a cluster using the spatial distribu- 
tion of the member galaxies. To determine the cluster shape, we employ the dispersion ellipse 



of the bivariate normal frequency function of position vectors (see, e.g.. iTrumpler &: Weaver 
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Carter & Metcalfe 


1980; 


Bureett et al. 


2004 


is defined bv Ixrumpler & Weaver 


Jl952 


) as 



times the maximum density of a set of points distributed normally with respect to two corre- 
lated variables, although the points need not be distributed normally in order to determine 
the proper cluster shape. From the first five moments of the spatial distribution, 



1 ^ 

N ^ " 



1 ^ 



i=l 
N 



N 



i=l 
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(22a) 
(22b) 
(22c) 
(22d) 



where X and Y are clustercentric distances in the direction of right ascension and decli- 
nation, respectively, the semimajor and semiminor axes of the ellipse, T a and F^, are derived 
by solving the equation 

F2 1^11 
/i02 - F 



A*20 



/ill 







(23) 



The position angle of the major axis, measured from north to east, is given by 



92 = cot"^ 
and the ellipticity is defined by 



/i02 



/^ll 



l-I^. 



+ 



(24) 



(25) 



The major and minor axes, position angles, and ellipticities derived for our sample 
clusters are listed in Table HJ and the dispersion ellipse is shown in Figure [31 The ellipticities 
are in the range e = 0.15 — 0.36, indicating no strong elongation. 
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5.2. Analysis of Substructure 

The analysis of substructure is a useful diagnostic tool for understanding the dynam- 



ical state of g a. 



axy clusters. A useful discussion of s e yeral substructure tests is given in 



Pinkney et al.l (Il996l ). As described in iHwang &: Led (120071 ). we derived number density 



maps using the spatial position (two-dimensional; 2D) of member galaxies, and performed 
one-dimensional (ID) and three-dimensional (3D) substructure tests for our sample clusters. 

The majority of ID (velocity histogram) substructure tests are n ormality tests. We 



present the results of five ID tests in Table [51 The values of the I-test (ITeague et al.lll990l ) 
are shown in the second and the third column; /go is the critical value for rejecting the 
Gaussian hypothesis at 90% confidence. Therefore, a velocity distribution is considered to 
be non-Gaussian if / > Iqq. We find that four of the sample clusters (A779, 1650, 1651, and 
2199) do not satisfy the Gaussian hypothesis using the I test. 

The skewness, which is a measure of the degree of asymmetry of a distribution, and the 
confidence level at which it rejects normality are given in the fourth and the fifth columns 
of Table [5l Positive or negative skewness indicates that the distribution is skewed to the 
right or left, respectively, with a longer tail to one side of the distribution maximum. The 
kurtosis, which is the degree to which a distribution is peaked, and the confidence level at 
which it rejects normality are given in the sixth and the seventh columns. Positive values 
indicate pointed or peaked distributions, while negative values indicate fiattened or non 
peaked distributions. The skewness test rejects a Gaussian distribution with a confidence of 
over 90% for A85, 1651, and 2199. The kurtosis test rejects the hypothesis of Gaussianity 
for A779, 1650, 1651, 2034, and 2199 with confidence of over 90%. 

From the eighth to the eleventh colu mn of Table [5l we p resent the asymmetry index 



(AI) and the tail index (TI) introduced by iBird fc Beerd (119931 ) along with their confidence 



levels. The AI measures the symmetry in a population by comparing gaps in the data on 
the left and right sides of the sample median, and TI compares the spread of the data at 
the 90% level with the spread at the 75% level. The Gaussian hypothesis for A779, 1800, 
and 2199 are rejected by the AI test, and that for A85, 1650, 2199 by the TI test with a 
confidence level of over 90%. 

We have constructed number density contours for the clusters using different bin sizes 
of 0.25R' X 0.25R' Mpc^ depending on the cluster size [R' = 4(rArB) 1/2/1.5]. and F^ are 
in units of Mpc, 4 is an arbitrary constant, and 1.5 (Mpc) is a normalization constant. We 
have smoothed the contours using a cubic convolution interpolation method. The contour 
interval, (max density in cluster)/6, is also determined according to the maximum number 
density of the clusters. We plot the number density contour map in the first and third 
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columns of Figure 



Usin g the velocity data and posi tional information on the galaxies, we have performed 
a A-test (jPressler &: ShectmanI Il988l ). which computes local deviations from the systemic 
velocity (fsys) and dispersion (dp) of the entire cluster. For each galaxy, the deviation is 
defined by 



where A^„„ is the number of galaxies that defines the local environment, taken to be 
~ ^ga/^^ in this study. The sum of b over all galaxies in a cluster. A, is used to quantify 
the presence of substructure. It is approximately equal to the total number of galaxies in a 
cluster in the case of no substructure, while it is larger than the total number of galaxies in 
a cluster in the presence of substructure. 

The statistical significance of the deviation is computed by Monte Carlo simulations. 
Velocities are randomly assigned to the galaxies at their observed positions, and Agim is 
computed for each simulated cluster. We construct 1000 simulated clusters and compute 
Agim for each simulation. We present Aobs which is computed using real data and the 
fraction of simulated clusters with Asim > Aobs in the final two columns of Table [51 Small 
values of /(Agim > Aobs) indicate statistically significant substructure. A779, 1795, 2199, 
and 2734 have much smaller values than other clusters, indicating significant substructures. 

We plot the positions of cluster galaxies, represented by circles with radii proportional 
to e"^, in the second and the fourth columns of Figure [161 A large circle denotes a galaxy that 
is deviant in either velocity or dispersion compared with nearby galaxies; therefore, groups 
of large circles indicate the presence of substructure. No strong substructures in A85, 1650, 
1651, 1800, 2034, and 2670 were found with the /(Asim > Aobs) test, as confirmed in Figure 
[T6l In the number density contour maps of Figure [161 we present the galaxies with b < 2.0 
that were used for orbit analysis (dots) and those with 6 > 2.0 in the substructure (crosses) 
separately, in order to show the usefulness of 6 parameter to identify the substructure. 



5.3. Dynamical Status of the Clusters 



A study of the dynamical state for our sample clusters is useful for interpreting the 
galaxy orbits in clusters. It is expected from self-similar models that a relationship oc CTp 
between the X-ray lumin osity and the velocity dispersi on of the clusters w ill exist (e.g. 
Quintana fc Melnicklll982l ). For galaxy clusters in SDSS, [Popesso et al.l (120051 ) found oc 
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^3.68±o.25 galaxy clusters in 2dFGRS. iHilton et all (120051 ) found a relation oc af, 



r4.8±0.7 



and suggested that high-Lx clusters are more dynamically evolved systems than the low-Lx 
clusters. In Figure [TTl we plot the X-ray luminosity as a function of velocity dispersion and 
virial mass for our sample clusters {filled circles) for which X-ray luminosities are available 
in the literature, compared with other clusters {open circles) out of the 113 selected clusters. 
We plot the X-ray luminosities from different literature sources separately. It shows that 
correlates well with Up and Myir. The power-law slopes are found to be in the range 3.7 — 4.9 
depending on X-ray references, which is consistent with the results of previous studies. Our 
sample clusters except for A779 are usually X-ray bright compared with the other clusters, 
because it is easy to derive X-ray mass profiles for the X-ray bright clusters. Interestingly, 
A779 and A1795 show a significant deviation from the relation L^-ap and Lx-Mvir. 

The dynamics of the brightest cluster galaxies (BCGs) or cD galaxies in galaxy clusters 



are a lso useful for understanding the formation history of the clusters (see, e.g.. lOegerle &: Hill 



20011 ). In particular, the peculiar velocity of the BCG, defined by Vp = vbcg — Vd, where 



vbcg is the observed velocity of BCGs and f ni is the mean velocity of the cluster, is a useful 



indicator of the dynamical state of a cluster (jOegerle fc Hillll200ll ). To estimate the pecu 



liar velocities of BCGs in clusters, we first identified the BCGs that have the brightest 6j 
magnitudes in the catalog of member galaxies for the 113 selected galaxy clusters. Then we 
conducted a visual inspection of cluster images to determine whether there are any galaxies 
brighter than the selected BCGs, using the catalog of member galaxies. Since some very 
bright galaxies in clusters were not covered as a result of observational difficulties such as 
fiber collision and saturation, we finally selected 24 galaxy clusters for which the BCGs in 
the catalog agree with those in the images. Using this sample of 24 galaxy clusters, we 
plot, in Figure [T8| the cluster velocity dispersion as a function of the absolute value of the 
BCG peculiar velocity, the absolute value of the BCG peculiar velocity as a function of 
clustercentric distance, and the redshift of the clusters as a function of the BCG absolute 
magnitude in the 6j band. 

It appears that the absolute values of the peculiar velocities of the BCGs are in the 
range 5—802 km s~^ and that the median value of the peculiar velocities is 165 km s~^. 
The mean uncertainty on the values of the peculiar velocities is 117 km s~^. In addition, 
the clustercentric distances of the BCGs are in the range 0—350 kpc, and the median value 
of the clustercentric distances is 75 kpc. Interestingly, the BCG peculiar velocities for two 
(A1651 and A2670) clusters are larger than that median value, indicating dynamical non- 
equilibrium. 
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6. Discussion 

We summarize the global kinematic properties for our sample clusters in Table [61 The 
first column gives the Abell identification. The second through fourth columns list the 
existence of substructure as indicated by the ID, 2D, and 3D tests, respectively : "Yes" 
is assigned if three of five tests indicate the presence of substructure for ID test, if the 
substructure at the contour level of 2x(max density in cluster)/6 is seen for 2D test, and if 
f{Asim > Aobs)< 0.1 for 3D test. The cluster morphology determined in §5. H is given in the 
fifth column: "spherical" for ellipticity less than 0.2 and "elongated" for ellipticity greater 
than 0.2. "?" is given for A85 since the spatial coverage is incomplete. The dynamical 
status determined in §5.31 is given in the sixth and seventh columns. For the scaling relation 
between the X-ray luminosity and the velocity dispersion or the virial mass, "No" is assigned 
if the deviation of X-ray luminosities from the scaling relations derived in this study is larger 
than one standard deviation in at least four out of eight panels in Figure [T71 For the peculiar 
velocities of the BCGs, "No" is assigned if the significance of the peculiar velocity of the 
BCG, S [= |fp|/(eBCG + ^ci)"'^''^]' is larger than 3. eBCG is the measurement error of the BCG 
velocity and ed is defined by ed = a^lNgai, where ad is the velocity dispersion of the cluster. 
The description on the result for the individual cluster is given Appendix A. 



6.1. Effects of Different X-ray Mass Profiles 



Since we used one X-ray mass profile per one cluster among several X-ray mass profiles, it 
is important to examine the results using different mass profiles. Moreover, the mass profile of 
Reiprich fc Bohringerl (120021 ) that we used in §3.11 adopted a constant gas temperature model 
to derive X-ray mass profile, but recent X-ra y observations showed that gas temperature is 
not constant over the clustercentric distance ( Markevitch et al.lll998l : iDe Grandi fc Molendi 



2OO2I : iPratt et al.l 120071 ). The isothermal model is known to make the mass profile steeper 



than that of non-isothermal mo dels, leading to underesti mate the cluster mass at small radii 
and overestimate at large radii. iMarkevitch et al.l (119981 ) showed that the isothermal model 
leads to underestimate the cluster mass by a factor of 0.74 of the non-isothermal model 
at one core radius and overestimate by a factor of 1.43 of the non-isothermal model at six 
core radius. To investigate the effects of different mass profiles on the determination of 
galaxy orbits, we compute /3orb in each cluster with various mass profiles, and present the 
results in Figure [T9l The ov erestimation of the cluster mass from the isothermal model of 
Reiprich fc Bohringerl (120021 ) is expected to give larger value of calculated VDP than the true 
value, and the value of [3r,rh would decrease . Indee d it is seen that the cluster masses of the 
isothermal model in iReiprich fc Bohringeii (120021 ) are larger than those of non-isothermal 
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model in ISanderson et al.l (120031 ) at r ~ r2oo for A1651 and A2199 (F i g- El) - Therefore, 
the calculated /Sorb with the isothermal mo del in Reiprich &: Bohringerl (120021 ) are smaller 
than those with the non-isothermal model ISanderson et al.l (120031 ) [Fig. [12] (b) and (c)]. 
However, the orbit determination does not significantly change, because these values of /5orb 
with various mass profiles are not much different and accord within the uncertainty. It is 
noted that the gal axy orbits can be det ermined differently as s een in the case of A21 99 
(isotropic orbit for ISanderson et al.ll2003l . but tangential orbit for iMarkevitch et al.lll999l at 
r ~ 0.35 Mpc; see Appendix lA.Sp . 

The effect of the most accurate mass profile obtained with Chandra is seen in (d), (e), and 
(f). Since the Chandra mass profile is available for only A1795 and the discrepancy among the 
mass profiles in A1795 is not large, the calculated /5orb with Chandra mass profile appears to 
accord with those from other mass profiles within the uncertainty. The discrepancy becomes 
larger at large radius (r > rx), but the reliability at this region is low because the mass 
profile is determined from the extrapolation of that for the inner region. 



6.2. Effects of Different Galaxy Samples 



Since significant fractions (40 — 70%) of the clusters show substructures, in dicating that 
they are in the process of merging (jjones &: Formanlll999l : iRamella et al.l 120071 ). it is impor- 
tant to study the effect of the presence of substructure on the determination of galaxy orbits. 
It is expected that the presence of substructure leads to an increase of the observed VDP 
and the change of the galaxy number density profile. As a result, the orbit determination 
can be diff erent depending on the de gree of inclusion of the galaxies in the substructure. 
Previously, iBiviano fc KatgertI (120041 ) showed that the galaxies in the substructure appear 
to have tangential orbit as a whole. Therefore, the orbit analysis including the galaxies in 
the substructure may make /5orb have low values. 

Figure [2D] shows the effect of the presence of substructure on the determination of galaxy 
orbits. As seen in (b) — (e), the calculated /3orb including the galaxies in the substructure 
{dotted lines) tends to be lower values than those excluding the galaxies in the subst r uctur e 
{solid and dashed lines), being compatible with the results of iBiviano Sz KatgertI (120041 ). 
However, A85 [Fig. [20] (a)] and late-type galaxies in A1795 [Fig. [20] (f)] do not show 
similar behavior to other cases, since the changes in observed VDPs and in galaxy number 
density profiles are not similar to other clusters. In conclusion, the orbit determination 
can be changed significantly [e.g.. Fig [20] (f)], depending on the degree and the location of 
subclustering, but is not much changed for our sample clusters. 
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6.3. Comparison with the Previous Studies 



For 10 Abell clusters, we have found that the orbits of galaxies are consistent with 
isotropic orbits in most clusters, although it is difficult to conclude strongly for some clus- 
ters due to the large errors (e.g., A779, A1650) and the variation as a function of cluster- 
centric distance (e.g., A1795, A1800, A1795, A2199) in calculated Porh- Isotropic galaxy 
orbits for majority of our sample clusters , are consistent with t ho se for composite clusters in 
the previous larg e cluster surveys (e.g., ICarlberg et aDll997al jbl: | van der Marel et al.l l2000l : 
Rines et al.l 120031). Existence of aniso tropic galaxy orbits is also reported in other studies. 
For A2218, iNatarajan fc KneibI (119961 ) found tha t the galaxy orbits ar e tangential at small 
radii {R < 400 kpc) and are radial at large radii. iBenatov et al.l (120061 ) analyzed the galaxy 
orbits for five clusters up to r ~ r200; and reported that A2199 and A496 are consistent with 
tangential orbits, while A2390 radial orbits. Two clusters (MS1358 and A576) appear to 
have radial orbits, but these clusters may not be in hydrostatic equilibrium. 

Numerical simulations sho wed that cluster ga laxies are likely to follow isotropic orbits 
as seen in many clusters (e.g., iGhigna et al.lll998l ). In detail, isotropic orbits are preferred 
in the inner region, and radial orbits in the outer region, indicating existence of accret 
ing galaxies from the clus t er outskirt ( Crone et al.lll994l: ICole fc LacevI Il996l : iTormen et al. 



1997; Thomas et al. 1998; Colin et al. 



2000 



Faltenbacher et al. 



20051 ). However, tangential 



galaxy orbits found in the observational data, are not commonly seen in sim ulation data 
(ITormen et al.l 119971 : iThomas et al.lll998l ). Interestingly, iBenatov et al.l (120061 ) reported an 
offset in orbital profile for their observed and simulated clusters in the sense that the simu- 
lated clusters show preferentially tangential orbits. 

For the orbital difference between early-type and late-type galaxies, the previous studies 
using the composite cluster, le d to the conclusion that early-types have quasi-isotropic orbits 



and late-types radial orbits (IBiviano et al.l 119971 : iMahdavi et al.l Il999l : iBiviano fc Katgert 



20041 ). Interestingly, the clusters in which orbital difference among subsamples were studied 
in this study, show no significant difference between them. In order to investigate the origin 
of difference between this study and the previous studies, we show VDPs and projected 
galaxy number density profiles for several composite clusters in Figure [211 We construct 
the composite clusters for three different samples: Sample A - using four clusters in which 
orbital properties for subsamples are studied (A779, 1650, 1795, and 2199), Sample B - using 
six clusters in which orbital properties among subsamples are not studied (A85, 1651, 1800, 
2034, 2670, and 2734), and Sample C - using 62 clusters without definite merging evidence, 
selected in §2.21 All samples are constructed using the galaxies outside the substructure 
{5 < 2.0). 



The VDP for the composite cluster is derived in the combined distribution of velocity 
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relative to the cluster mean and normalized by the velocity dispersion of each cluster. It is 
known that the shape of VDP provides some hints for the galaxy orbits iii clusters, and varies 



depending on cluste rs (Iden Hartog fc Katgertlll996l : iMahdavi et al.lll999l : iMuriel et al.ll2002 



Aguerri et al.l 120071 ): flat, decaying, or rising VDP. It appears that the differences in VDPs 



and in projected galaxy number density profiles between early-type and late-type galaxies 
for Sample A are not significant at 0.3 < R/r2oo < 0.8 where the VDPs are available, while 
those in Sample B are large. In Sample C, the difference is significant with small errors: the 
values of velocity dispersion for late-types are larger than those for early-types, and the VDP 
for late-types is decreasing along the radius while that for early-types is nearly fiat. Sample 



C is c omparable to the composite cluster based on ENACS data used in iBiviano &: Katgert 



( 20041). and our resu lts for Sample C are similar to theirs. Thus, if we adopt the result of 



Katgert et al.l (120041 ) that early-types have isotropic orbits, then the orbits for late-types in 



our composite cluster (Sample C) can be regarded as radial orbits, which is consistent with 
the results of the previous studies. Similarly, since not only VDPs but also projected galaxy 
number density profiles are not distinguishable in Sample A at 0.3 < R/r2oo < 0.8, it is 
expected that the orbital difference between early-type and late-type galaxies in Sample A is 
not significant. If we use the galaxies with lower probability of belonging to the substructure 
(e.g., 6 < 1.4), the results change little, indicating the effects of the galaxies substructure is 
not significant. Therefore, no orbital difference between early-type and late- type galaxies in 
Sample A c a n be interpreted as their own characteristics of individual clusters. Previously, 



Rines et al.l (120051 ) also reported that the VDP of ELG is not different from the non-ELG 



for CAIRNS c l usters , but their spatial distributions are different. In numerical simulations. 



Diaferio et al.l (l200ll ) found that there is no significant difference in orbital properties between 



red and blue galaxies : both red and blue galaxies are close to being isotropic. Thus, it is 
necessary to study what makes the diversity of galaxy orbits using more observational and 
simulation results. 



7. Summary 

We present the results of a study for the galaxy orbits in galaxy clusters using a spec- 
troscopic sample of galaxies in SDSS and 2dFGRS. Our results are summarized as follows: 

1. We have determined the member galaxies of 731 and 230 Abell clusters covered in 
SDSS and 2dFGRS, respectively. We have selected 10 clusters using three criteria : 
the number of member galaxies is greater than or equal to 40, the spatial coverage is 
complete, and X-ray mass profile is available in the literature. 

2. For the selected 10 clusters, we derived the radial profile of the galaxy number density 
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and velocity dispersion for all, early-type, and late-type galaxies outside the substruc- 
ture [6 < 2.0). 

3. We have investigated the galaxy orbits for our sample clusters with constant and 
variable PoA through the clustercentric radius using Jeans equation. The resulting 
galaxy orbits based on two methods appear to be consistent. 

4. Using all member galaxies, most of our sample clusters are found to be have isotropic 
orbits, although it is difficult to conclude strongly for some clusters due to the large 
errors (e.g., A779, A1650) and the variation as a function of clustercentric distance 
(e.g., A1795, A1800, A1795, A2199) in calculated povh- For four clusters (A779, 1650, 
1795, and 2199) the orbital difference between early-type and late-type galaxies appears 
not to be significant. 

5. We have determined the cluster morphology using the dispersion-ellipse method, and 
have found the ellipticity in the range of 0.15 — 0.36, indicating no strong elongation 
for our sample clusters. 

6. We have investigated dynamical status for our sample clusters using substructure (ID, 
2D, and 3D) tests, the relation between X-ray luminosity and the velocity dispersion 
or the virial mass, and peculiar velocity of BCGs. It is found that the majority of our 
sample clusters are dynamically relaxed system. 
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A. Properties of Individual Clusters 
A.l. Abell 85 



A85 was covered in SDSS, but the outer region {R > 26') was not surveyed (see Fig. [3]). 
This cluster is in the complex of clusters Abell 85/87/89, and the member galaxies selected 
in this st udy include only Abell 85/87. The structure of this cluster is well studied by severa l 
authors Jhima Neto et al.ll200ll : lDurret et aDll998l : lKempner et allbooi burret et al.ll2005l l: 
a southern blob (or group) that is merging from south (~ 10' from the cluster center), a 
subcluster that is merging from southwest (~ 4' from the cluster center), an extended 4 Mpc 
X-ray filament that is probably made of groups falling onto the main cluster. Interestingly, 



X-ray studies base d on recent Chandra and XMM- Newton data (e.g., iKempner et al.l 12002 
Durret et al.ll2005l ). provided an evidence of past and present merger activity for this cluster. 
Our galaxy data showed that only A87 (~ 30', or 1.9 Mpc from the southwest of the main 
cluster) is detected as a substructure, but some of the galaxies in A87 with 6 < 2.0 were 
removed for the orbit analysis. In addition, a small deviation from the scaling relation 
between X-ray and optical data and a small peculiar velocity of the BCG do not imply 
dynamical non-equilibrium. The galaxy orbits are consistent with isotropic orbits within the 
error at 0.3 < r < 1.5 Mpc. Therefore, it appears that main cluster of A85 is stabilized at 
present, although it experienced a merging in the past, and is not disturbed significantly by 
the current merging event. We also find that the determined galaxy orbits do not change 
significantly, if we use different X-ray mass profiles as seen in Figure [13 
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A. 2. Abell 779 



A779 is one of the nearest (z~0.023) clusters and shows the faintest X-ray luminosity in 
our sample. The velocity dispersion is the smallest (cTp = 491^36 km s~^) among our sample, 
and the number of late- types {Ngai = 78) is larger than that of early- types {Ngai = 67). 
This cluster has a cD galaxy, NGC 2832, which is known to be at rest in cluster potential 
(lOegerle fc Hillll200ll ). but was not included in the spectroscopic sample of SDSS. The VDP 
using all member galaxies with 6 < 2.0 declines at R 5, 0-8 Mpc (see Fig. [9l), which is 
consi stent with the results of the previous studies (IMahdavi et al.l Il999l : iRines fc Diaferio 
20061 ). The galaxy orbits using all members appear to be consistent with radial or isotropic 
orbit within the uncertainty up to r ~ r2007 but it is difficult to conclude strongly due to 
the large error of calculated PoA- The orbits of early- types appear to be consistent with 
being radial or isotropic within the uncertainty through the radius, but strong conclusion is 
difficult due to the large error of calculated Porh- The orbits of late-types are consistent with 
being radial in the inner region (r < 0.1 Mpc) where the observed VDPs are not available, 
but with being isotropic within the error at r < 1 Mpc and become tangential at outer 
region. 



A.3. Abell 1650 



A1650 was included in both SDSS and 2dFGRS. A recent XMM-Newton observation 
showed no evidence for spat ial temperature variation and surface brightness irregularity, 
indicating a relaxed system (ITakahashi &: Yamashitall2003l ). Similarly, we found no strong 
evidence for non-equilibrium (see Table [6]). Figure [8] shows that Bj — Rp color of the BCG 
in this cluster is ~ 0.8 redder than the color expected from the OMR. However, this color 
deviation of the BCG in th is cluster is not seen in other studies (IMargoniner fc de Carvalho 
2000l : IPimbblet et al.l[2006l ). Therefore, the Bj - Rp color of the BCG given in 2dFGRS needs 
further investigation. The galaxy orbits for all member galaxies are found to be isotropic 
within the uncertainty or tangential at r > 0.1 Mpc. Interestingly, the orbits for early-type 
and late-type galaxies appear to be tangential at r > 0.3 Mpc. 



A.4. Abell 1651 



A1651 is included in both SDSS and 2dFGRS. This cluster is known to be a dynamically 
relaxed system: re gular shape in the ROS A T PSPC image and symmetric temperature profile 
from ASCA data (IMarkevitch et al.lll998l ). Interestingly, a large peculiar velocity {vp = 350 
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km s~^) of the BCG is found i n this study (see Fig. [T8l) . indicating a dynamically unrelaxed 
system. lOegerle fc Hilll (I2OOII ) already identified a peculiar velocity of the BCG (252 — 258 
km s~^), but they concluded that the peculiar velocity is not large since a significance, 5* is 
less than 3. They obtained 5* ~ 1.3 with 39 member galaxies. However, we obtain 5* = 3.46 
with a large number of member galaxies {Ngai = 258) and small velocity error, indicating a 
significant {S > 3) measured peculiar velocity. If we use the galaxies outside the substructure 
with 6 < 2.0, the result change little. This cluster was included in Las Campanas/ Anglo- 
Australian Telescope Rich Cluster Survey (LARCS: iPimbblet et al.ll2006l ). where a recession 
velocity of the cluster and the cD galaxy was determined to be 25466 ± 57 and 25571 ± 27 
km s~^, respectively. However, we obtain the mean velocity of the cluster (biweight location 
used in this study), Vd = 25278 ±95 km s~^ using their data, confirming significant peculiar 
velocity, Vp = 293 km s^^. In contrast, the galaxy orbits are found to be isotropic within 
the uncertainty through the radius (see Fig. [TTl) . which means that the galaxies are in 
equilibrium with the cluster potential, being consistent with the result from X-ray data. 
Therefore, it would be interesting to investigate what makes the large peculiar velocity of 
the BCG. 



A. 5. Abell 1795 



A1795 is the most X-ray-luminous cluster in our sample. A cD galaxy in this cluster 
was not included in the spectrosco pic sample of SP SS. A larg e pecuhar velo c ity (365 krn 



s-^) of the cD galaxy was found bv iHill et all fll988h . but later lOeeerle fc Hilll fll994l . l200lh 



reported the peculiar velocity is not large (150 — 180 km s ^) and its significance is less than 
3. Al though detailed Chandra data showed that the central core of t his cluster is not relaxed 



Markevitch et al.ll200ll : Ettori et all I2OO2I : iFabian et all 120011 1. o ther wide-field X-ray 



[e.g. 

data indicate that the cluster is close to being dynamically relaxed (IBriel fc Henryl Il996 



Buote fc Tsailll996l : iTamura et allboOlh . We found that the galaxy orbits derived from all 
member galaxies with S < 2.0 are consistent with being tangential at 0.3 < r < 1.2 Mpc, 
but being isotropic within the uncertainty in the outer region. Interestingly, the orbits of 
early-types appear to be tangential at r > 0.3 Mpc, while those of late-types to be tangential 
or isotropic within the error. The results with different mass profiles are shown in Figure! 



A.6. Abell 1800 



A1800 is one of the member clusters in the Bootes supercluster (lEinasto et al.ll200ll ). and 
has a relatively small number ratio (~ 0.16) of late- types to early-types among our sample 
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clusters. There is no strong evidence of dynamically non-equilibrium from the substructure 
tests and the scaling relation between X-ray and optical data. The galaxy orbits change from 
radial (r < 0.2 Mpc), to isotropic (0.2 < r < 0.8 Mpc) within the error and to tangential 
(0.8 < r < 1.8) Mpc. 



A. 7. Abell 2034 

A2034 is the most distant cluster (z~ 0.113), and has the largest velocity dispersion 
{UOSlljG km s"^) in our sample clusters. Interestingly, there are only four late-types out 
of 78 member galaxies. Since the observed VDP is larger than the upper envelope (/3orb = 
0.99) of the computed one over the radius, it is difficult to conclude using top and middle 
panels. However, the bottom panel indicates that the galaxy orbits are radial in the outer 
region (r > 1 Mpc), and those are consistent with isotropic orbits within the uncertainty at 
0.2 < r < 1 Mpc. 

Larger values of observed VDP than the strong radial orbit can be interpreted as non- 
equilibrium of this cluster, although substructure tests and the scaling relation between X-ray 
and optical data in this study imply no strong substructure. However, other studies showed 
direct evidence of dynamical non-equilibrium. For example, detaile d analysis of Chandra 



data revealed evidence for an ongoing merger (iKempner et al.ll2003l ): northern cold front 
that is a discontinuity of the surface brightness on the northeast edge of the cluster, large 
concentration of galaxies including a cD galaxy just ahead of the cold front, and excess of 
emission to the south of cluster. In addition, a cD galaxy in the main cluster is offset ~ 1' 
from the X-ray centroid. Unfortunately, the cD galaxy is not included in the spectroscopic 
sa mple of SPSS. If w' e adopt a receding velocity {vp = 33445 ± 35 km s~^) of the cD galaxy 



m 



Miller et al.l ((2002), then the peculiar velocity for this galaxy is 477 ± 163 km s~^ with 
the significance of 5* = 3.378. This large peculiar velocity becomes larger if we use recession 
velocity {vd = 34373 km s~^) of the cluster determined in the same reference. 



A.8. Abell 2199 



The number of member galaxies {Ngai = 754) in A2199 is the largest among our sample 
clusters, and the number {Ngai = 372) of early-types is comparable to t hat (Nggi = 382) 



of lat e-types. A2199 is known to be one of the richest, r egular cluster (jMarkevitch et al. 



1999f) . and contains several infalling bound subclusters (iRines et al.l l200ll ). iRines et al. 
2OO2I ) derived the VDP up to 8 Mpc using ~ 300 member galaxies in A2199, and found 
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that th e observed p r ofile i s consistent with an isotropic orbit based on the caustic mass 
profile. iLokas et al.l (120061 ) . using 180 member galaxies within i? ~ 1.1 Mpc, determined 
Porh = ~0.55t2:75 based on the joint analysis of velocity dispersion and kurtosis, which is 
consistent with a tangential or isotropic orbit. We found, using an independent X-ray mass 
profile and a smoothed velocity profile, that the galaxy orbits are isot ropic at 0.1 < r < 0.7 
Mpc within the error, and are radial in the outer region. However, iBenatov et al.l (120061 ) 
found that A2199 has ta ngentially anisotropic orb it in the inner region {R < 1.1 Mpc), using 
an X-ray mass profile of iMarkevitch et al.l (119991 ) and polynomial fit of V P P. This discrep 
ancy i s demonstrated in Figure [T9l (c), and is because the mass profi le of iMarkevitch et al. 



( 119991 ) they used is larger than that we used (ISanderson et al.l l2003l ) (see ^ for more de- 
tail). We also investigated the orbital difference between early-type and late-type galaxies in 
A2199. Interestingly, both subsamples appear to have radial through the radius, indicating 
they infall from the outer region. 



A.9. Abell 2670 



Abell 2670 was included in both SDSS and 2dFGRS. We found no strong evidence of 
substructure from the substr uctur e tests, but identified a large peculiar velocity of the BCG, 
which was already found by iBirdI (Il994a[). Although X-ray morpho l ogy see n by RO SAT is 
regular, indicating a relaxed system (IHobbs fc Willmorelll997l ). iBirdI (Il994bl ) concluded that 
this cluster may c onsist of four subclus ters that are merging along the line of sight, using 
~ 230 galaxies in ISharples et al.l (119881 ) catalog. iHobbs fc Willmord (119971 ) also suggested 
a merging activity since the observed VDP is much greater than that expected from the 
X-ray mass profile. In contrast, we found isotropic galaxy orbits, which imply a relaxed 
system. The observed dispersio n profile in this study is not significantly different from 
that in lHobbs fc Willmord (Il997l). However , X-ra y mass profiles adopted in this study are 



larger than that us ed in 



Hobbs&Wil 



niord (119971) . thus the value of expected dispersion is 
larger than that in iHobbs fc Willmord (119971 ). Therefore, other evidence (e.g., irregular gas 



temperature map) is needed to confirm the merging activity of this cluster. 



A. 10. Abell 2734 



Burgett et al.l (120041 ) found no strong substructure using 1 25 velocity da t a from the 



2dFGRS data . Sim ilar results were found by other authors: ISolanes et al.l (Il999l ) and 



Biviano et al.l (120021 ) identified no substructure using 45 and 77 member galaxies in ENACS 
data, respectively. However, our 3D substructure test resulted in a strong substructure in 
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this cluster. It is due to the larger spatial coverage than previous studies, therefore, the sub- 
structure at i? ~ 30' (~ 2 Mpc) to the south-east is included in this st udy. Comparison of 



ROS A T X-reiy image with the optic al one also shows a good agreement (IKolokotronis et al. 



200ll ). Iden Hartog fc Katgertl (119961 ) derived a flat VDP up to R~ 1.5 Mpc using 77 velocity 



data in EN ACS. The galaxy orbits appear to be isotropic within the uncertainty through the 
radius, indicating a dynamical equilibrium, which is consistent with the previous studies. 
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Table 1. The Sample of Galaxy Clusters 



All Early-types Latc-typcs 

R.A. Decl. B-M cz o> A^gai o> Ng^i o> 'W~i 

Cluster (J2000) (J2000) Type Survey* (km s"!) kpc/arcmin (km s"!) (km s"!) (km s"!) 



A0085 


00 


41 50.09 


-09 18 06.8 


I 


S 


16545l|i 


64.3 


926+51 


208 


804+^1 


170 


1114ll»« 


38 


A0779 


09 


19 41.28 


+33 45 46.8 


I-II 


S,2 


6947l|i 


28.1 




145 


471+73 


67 


479115 


78 


A1650 


12 


58 41.09 


-01 45 24.8 


I-II 


8,2 


251871^^ 


94.7 


762+39 


258 


7221^1 


168 


829l^f, 


89 


A1651 


12 


59 21.50 


-04 11 40.9 


I-II 


S 


25222^^^ 


94.8 


88lll« 


258 


828+^1 


169 


996 7g 


84 


A1795 


13 


48 52.97 


+26 35 44.2 


I 


S 


18775l«^ 


72.4 


755l« 


163 


729+*^ 


127 


806±™ 


36 


A1800 


13 


49 21.65 


+28 06 13.0 


II 


S 


22660l|5 


86.0 


7621S 


93 


722l« 


80 




13 


A2034 


15 


10 11.74 


+33 30 52.9 


II-III 


s 


33922tjf° 


123.3 


1208t7g 


78 


11851^^ 


74 


174611^2^ 


4 


A2199 


16 


28 37.97 


+39 32 55.3 


I 


s 


922ltg 


36.9 


726122 


754 


703+29 


372 


75015J 


382 


A2670 


23 


54 13.39 


-10 24 46.1 


I-II 


8,2 


22837lg2 


86.6 


8441^^ 


106 


840lt^ 


82 


870li;j^ 


24 


A2734 


00 


11 20.71 


-28 51 18.0 


III 


2 


18256tfg 


70.5 


lOOllli 


192 


83lt^° 


120 


1175t?2 


68 



* 'S' is covered by 8DS8 and '2' by 2dFGRS. 
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Table 2: X-ray Properties for Our Sample Clusters 





Lx(0. 1-2.4 keV)'' 





' c 




T(0) 


'y 




Cluster 


/iq44 QT-rrc prn"^ S^^") 












Ref.'' 


Ann85 


Q 789 


76 


293 8 




10 97 


1.32 


S03 


Au My 


u.uyu 


U.o4 






d.DO 






A1650 


7.308 


0.78 


214.2 




9.73 


1.19 


S03 


A1651 


8.000 


0.70 


192.5 




6.87 


1.10 


S03 


A1795 


10.124 


0.83 


288.4 




9.98 


1.17 


S03 


A1800 


2.840 


0.77 


280.0 


4.02 






R02 


A2034 


6.850 












D03 


A2199 


4.165 


0.60 


76.9 




4.17 


1.15 


S03 


A2670 


1.469 


0.55 


83.8 




5.99 


1.04 


S03 


A2734 


2.365 


0.62 


151.4 


3.85 






R02 



^ A0779 flEbeling et al.ll2000l). A2034 f|Ebeling et al.lll998r) . A2670 fiBohringer et al.ll2004j ). 
& othe r clusters (Reiprich fc BoliringeiT2002 '). 

^ R0 2 jReiprich fc Bohringei]l2002h . S03 f|sanderson et al.ll2003h . and D03 JPemarco et al. 
2003h . 



Table 3. Galaxy Number Density Profile for the Clusters 





I'O.NFW 


I's.NPW 


I'O.Hcr 


f's,Hcr 


i^O.KEKOS 


''s,KEK05 


''yg.KEKOS 


Cluster 


(arcmin 




(arcmin) 


(arcmin^'^) 


(arcmin) 


(arcmin 




(arcmin) 




All Galaxies 


A0085 


1.34£; 


- 3 


30.2 


6.85E - 4 


54.6 


1.79E 


- 2 


15.6 


9.0 


A0779 


8.82£; 


- 3 


7.4 


8.01E - 4 


29.1 


2.05B 


- 2 


16.0 


12.7 


A1650 


4.33_B 


- 3 


15.3 


1.59£; - 3 


33.6 


4.14E 


- 2 


13.0 


10.7 


A1651 


3.95E 


- 3 


15.3 


1.86E - 3 


30.1 


4.67E 


- 2 


10.4 


10.1 


A1795 


3.88E 


- 3 


13.7 


imE - 3 


28.4 


2.83E 


- 2 


16.4 


11.8 


A1800 


7.88E 


- 3 


8.0 


2.99S-3 


17.6 


5.96E 


- 2 


9.0 


11.5 


A2034 


1.88E 


- 2 


6.3 


6.89B-3 


14.0 


1A7E 


- 1 


6.3 


11.1 


A2199 


5.33E 


- 4 


40.5 


1.56S-4 


98.1 


8.07E 


- 3 


23.3 


9.9 


A2670 


3A3E 


- 3 


15.6 


1.61B- 3 


29.7 


2.mE 


- 2 


10.9 


9.1 


A2734 


2.35E 


- 2 


5.6 


6.81B- 3 


14.4 


1.52E 


- 1 


6.6 


11.7 


Early-type Galaxies 


A0779 


7.73E 


- 2 


2.3 


1.26E - 3 


17.3 


1.22E 


- 2 


16.1 


13.6 


A1650 


9.18E 


- 3 


9.1 


2.56E - 3 


22.9 


3.91E 


- 2 


13.2 


11.7 


A1795 


3.03E 


- 2 


5.2 


5.66E - 3 


15.6 


iA7E 


- 2 


15.4 


13.4 


A2199 


5.60E 


- 2 


4.7 


2.37E - 3 


24.5 


3.29E 


- 2 


18.6 


13.2 


Late-type Galaxies 


A0779 


8.96E 


- 6 


124.4 


5.29E - 6 


208.0 


2A1E 


-4 


28.8 


7.4 


A1650 


i.38E 


- 5 


119.6 


2.54B - 5 


200.1 


2.09E 


- 3 


17.0 


6.9 


A1795 


2.16E 


- 3 


9.9 


2.90E - 4 


32.9 


7.82E 


- 3 


14.8 


11.4 


A2199 


4.12E 


- 5 


105.8 


2.19E - 5 


189.8 


1.20E 


- 3 


26.7 


8.1 
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Table 4: Morphological Parameters of the Clusters 



Cluster 


Ta 




02 


e 




(kpc) 


(kpc) 


(deg) 




A0085 


135811 


868lt? 


1541^ 


'-'•'J"-0.05 


A0779 


11691^^ 




1611^ 


n 29+°-°^ 

'-'•^^-0.07 


A1650 




14201^^ 


12311 


n S4+0-05 

U-<J'±-0.05 


A1651 


211913^^ 


17981^^ 


52+^1 

'-'^-12 


n 1 5+°-°^ 

'^•-'-^-0.05 


A1795 


13651^^ 


11511^? 


47+14 

^' -16 


A 1 f^+0.06 


A1800 


10841^1 


910li 


63+}? 


•J- -'-'-'-0.06 


A2034 


123711 


9021^^ 


231}° 


rv 97+0.09 
-0.08 


A2199 


2182l-^;^ 


16461^2 


20+3 


'-'•^'-'-0.03 


A2670 


74411} 


5911^? 


401}^ 


'-'•^-'--0.06 


A2734 


14131^^ 


11521^1 


841^^ 


•J- -"-0-0.05 



Tab le 5: Parameters for Testing the Presence of Substructure 



Cluster I Iqq Skewness clRej Kurtosis clRej AI clRej TI clRej ^obs K^sim > ^obs) 



A008B 


1.03 


1.03 


-0.30 


93.0 


0.38 


81.4 


-0.22 


31.0 


1.18 


96.8 


251 


0.304 


A0779 


1.09 


1.04 


0.08 


35.6 


0.84 


95.2 


-1.01 


95.6 


1.16 


89.2 


245 


0.000 


A1650 


1.04 


1.02 


0.15 


71.4 


0.53 


91.2 


-0.43 


58.6 


1.21 


99.2 


326 


().15(i 


A1651 


1.05 


1.02 


-0.31 


96.4 


0.66 


95.2 


-0.61 


77.4 


0.99 


7.2 


345 


0.075 


A179B 


0.96 


1.03 


-0.04 


16.0 


-0.34 


51.4 


-0.17 


24.4 


1.03 


28.2 


246 


0.001 


A1800 


0.97 


1.05 


0.29 


77.2 


-0.54 


66.6 


0.87 


90.6 


0.87 


84.0 


118 


0.034 


A2034 


0.95 


1.06 


0.18 


53.0 


-0.82 


93.0 


0.34 


52.8 


0.95 


29.2 


94 


0.059 


A2199 


1.04 


1.01 


0.44 


99.8 


0.37 


92.2 


1.45 


99.2 


0.93 


92.6 


1456 


0.000 


A2670 


0.94 


1.05 


-0.01 


3.4 


-0.66 


86.0 


-0.15 


21.0 


0.94 


46.6 


102 


0.729 


A2734 


0.97 


1.03 


-0.07 


29.4 


-0.37 


65.8 


0.39 


58.8 


1.07 


62.2 


341 


0.000 



Table 6. Summary of Global Kinematic Properties for the Clusters 



Substructure? Equilibrium? 



Cluster 


ID 


2D 


3D 


Morphology 


X-Ray 


Vpec 


A0085 


No 


Yes 


No 


Elongated? 


Yes 


Yes 


A0779 


Yes 


No 


Yes 


Elongated 


No 




A1650 


Yes 


Yes 


No 


Elongated 


Yes 


Yes 


A1651 


Yes 


Yes 


Yes 


Spherical 


Yes 


No 


A1795 


No 


No 


Yes 


Spherical 


No 




A1800 


No 


Yes 


Yes 


Spherical 


Yes 




A2034 


No 


No 


Yes 


Elongated 


Yes 




A2199 


Yes 


Yes 


Yes 


Elongated 


Yes 




A2670 


No 


Yes 


No 


Elongated 


Yes 


No 


A2734 


No 


Yes 


Yes 


Spherical 


Yes 
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Fig. 1. — (a) Comparison of eclass parameter in SDSS with rj parameter in 2dFGRS, (b) 
histogram for rj parameter, and (c) histogram for eclass parameter. Sohd hne indicates 
the best hnear fit for eclass (< 0.09) and r] (< 1.1). The dotted vertical hne denotes 
the division value rj = —1.4 in 2dFGRS for galaxy subsamples, and dotted horizontal line 
indicates eclass=— 0.0640 in SDSS, which is equivalent to the division value rj — —1.4 in 
2dFGRS. 
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Fig. 2. — Radial velocity vs. clustercentric distance of galaxies and the velocity distribution 
for our sample clusters. Filled circles indicate the galaxies selected as cluster members, 
while open circles the galaxies not selected as cluster members. The horizontal dot-dashed 
lines indicate the systemic velocity of the clusters determined in Table [H The velocity 
distributions for the member galaxies are shown by hatched histograms, and those for all of 
the observed galaxies by open histograms. 
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Fig. 3. — Spatial distribution of cluster galaxies with measured velocities for our sample 
clusters. Early-type galaxies are plotted with crosses, while late-type galaxies are with open 
circles. All the observed galaxies are represented by dots. The ellipses indicate the (twice 
enlarged) dispersion ellipse, and solid lines denote the major and minor axis of the dispersion 
ellipse. 
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Fig. 4. — Distribution of the substructure parameter 6 for our sample clusters. The dotted 
vertical lines indicate the selection criteria for the galaxies in the cluster main body. 
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Fig. 5. — Total mass p rofiles for our sample cl usters. Dot-d a shed lines from 



Reiprich fc Bo 



iringerl (20021. R02), solid lines are from ISanderson et al.l (2003, S03 ), dot- 



ted lines from iDemarco et al.l ( 20031. DOS), dashed line s from IVikhlinin et al.l (120061, V06), 
and dot-dot-dashed lines from Markevitch et al. ( 19991, M99). The vertical dot-da shed and 
solid lines indicate the outer significance radius, r^, in lReiprich &: Bohringer J2002h and r2oo 
radius computed in this study, respectively. 
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Fig. 6. — Projected number density profiles of all (top panel), early-type (middle panel), 
and late-type (bottom panel) galaxies with 6 < 2.0 for A85, 779, 1650, 1651, 1795, and 1800. 
The solid, dashed, and dot-dashed lines indicate the projection of the best fits using NFW, 
Hernquist, and KEK05 profiles, respectively, for each sample. 



- 46 - 



1.00 




All 


0.10 






0.01 


A2034 




1.00 




Early 


0.10 


X 


X 


0.01 






1.00 




Late 


0.10 






0.01 








1 


10 


1.00 




All 


0.10 






0.01 


A2670 




1.00 




Early 


0.10 


X 


xx 


0.01 




X 


1.00 




Late 


0.10 






0.01 








1 


10 



01 



0.100 
0.010 
0.001 

0.100 
0.010 
0.001 

0.100 
0.010 
0.001 





All 


A2199 






Early 








Late 







10 



100 



1.000 




All 


0.100 
0.010 
0.001 


A2734 




1.000 
0.100 


7 X 


Early 


0.010 
0.001 




X 


1.000 




Late 


0.100 






0.010 
0.001 




® ® « 



R [arcmin] 



10 

R [arcmin] 



100 



Fig. 7.— Same as Fig. El but for A2034, 2199, 2670, and 2734. 
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Fig. 8. — Color-magnitude diagrams for our sample clusters. Filled circles indicate the 
photometric sample of galaxies within the i?max/2 from the cluster center, while open circles 
the spectroscopically selected member galaxies in the same region. Solid lines indicate the 
best linear fit of the CMR for early-type galaxies. BCGs are marked by open squares. 
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Fig. 9. — Velocity dispersion profiles of all (top panel), early-type (middle panel), and late- 
type (bottom panel) galaxies with S < 2.0 for A85, 779, 1650, 1651, 1795, and 1800. Filled 
circles indicate the velocity dispersion (up) at each point. The dispersion is calculated using 
the galaxies within the moving radial bin that is represented by a horizontal errorbar in 
each panel. The dotted lines denote 68% confidence interval on the calculation of velocity 
dispersion. The dashed horizontal line indicates the global value of velocity dispersion using 
all galaxies in each cluster. 
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Fig. 10.— Same as FigU but for A2034, 2199, 2670, and 2734. 
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Fig. 11. — Velocity dispersion profile (top panel) and the aperture VDP (middle panel), 
velocity anisotropy profile (bottom panel) for all galaxies of A85, 779, 1650, and 1651. 
Filled circles in top panels represent the measured VDP shown in Fig. [9], and those in 
middle panel denote the measured aperture VDP. Associated, dotted lines represent 68% 
confidence interval on the calculation of velocity dispersion. Three smoothly curved lines in 
top and middle panels, from a radial anisotropy to the tangential one (from top to bottom, 
Porh= 0.99, 0, and —99), represent the calculated dispersion profiles using galaxy number 
density profile of NFW (solid lines), Hernquist (dashed lines), and KEK05 (dot-dashed lines). 
Velocity anisotropy profiles (VAPs) determined in §4.21 in the bottom panel are represented 
by solid (NFW), dashed (Herquinst), and dot-dashed line (KEK05). The errors of VAPs 
for NFW are shown by dotted l ines. The vertical dot - dashe d and solid lines indicate the 
outer significance radius, rx, in iReiprich &: Bohringerl (120021 ) and r2oo radius computed in 
this study, respectively. 
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Fig. 12.— Same as Fig. [IH but for A1795, 1800, 2034, and 2199. 
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Fig. 13.— 



Same as Fig. [IH but for A2670 and A2734. 
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Fig. 14.— 



Same as Fig. [TT], but for early-type galaxies in A779, 1650, 1795, and 2199. 



-54- 



6 

o 
b 



800 

600 

400 

200 


800 

600 

400 

200 








A0779- 

'~-r~?.- • • ■ |^ 










_ ■ * f ">" V V * 


LVj r'V' i*; 







0.5 1.0 1.5 
R [Mpc] 





1 























.D 
O 


-1 












-2 
-3 








/a 




0.0 


0.5 


1.0 


1.5 



6 

a 
b 



1000 



1000 




0.0 

1 





Al?50" 


• — ■ 


•';■.■■.!•',■.■; 




~i — ^ 




' 1 

i 


_ . .. t 


■■■'""\ 




i " 

i 



0.5 



1.0 1.5 2.0 
R [Mpc] 



r [Mpc] 





0.0 0.5 



1.0 1.5 2.0 
R [Mpc] 



2.5 





1 ( 




- 






.D 

u 

O 


-1 : 




-2 




-3 : 









71 



0.0 0.5 



1.0 1.5 
r [Mpc] 



2.0 2.5 



6 

a 
b 



1500 

1000 

500 


1500 

1000 

500 







A2199" 






- i 









0.5 1.0 1.5 2.0 2.5 
R [Mpc] 




0.0 0.5 



1.0 1.5 
r [Mpc] 



2.0 2.5 



Fig. 15. — Same as Fig. [TT], but for late-type galaxies in A779, 1650, 1795, and 2199. 
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Fig. 16. — Left panels: Galaxy number density contour maps for our sample clusters. The 
member galaxies with 6 < 2.0 and 5 > 2.0 are represented by dots and crosses, respectively. 
The number density contours are overlaid. The plus signs indicate the cluster centroids, 
and the thick horizontal bars represent a physical extent of 1 Mpc. Right panels: Dressler- 
Shectman plots for the same clusters. Each galaxy is plotted by a circle with diameter 
proportional to e^. North is up, and east is to the left. 
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Fig. 17. — X-ray luminosity (0.1—2.4 keV) vs. the velocity dispersion and virial mass for 
the our sample clusters (filled circles) compared with the other clusters (open c i rcles) for the 
113 selected clusters. X-ray lumino s ities in (a) a nd (b) are from lBohringer et al.l (120001). those 



i n (c) and (d) from lEbeling et al.l (119981 . |2000| ). thos e in (e) and ( f ) frorn iBohringer et al. 



(120041 ). and those (0.5-2.0 keV) in (g) and (h) from iLedlow et all (l2003f ). The solid lines 
indicate the best fit for each panel. 
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Fig. 18. — (a) Velocity dispersion as a function of the absolute value of the pecuhar velocities 
of the BCGs, (b) absolute value of peculiar velocity as a function of the clustercentric dis- 
tances of the BCGs, and (c) redshift as a function of the absolute magnitudes of the BCGs 
in the 6j band for our sample clusters (filled circles) compared with the other clusters (open 
circles) among the 113 selected galaxy clusters. Histograms of the peculiar velocities (d), 
the clustercentric distances (e), and the absolute magnitudes in the 6j band (f) of the BCGs 
are also shown for the clusters including our sample clusters. The clusters in our sample are 
marked by filled circles in (d), (e), and (f). 
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Fig. 19.— VAPs for all galaxies of A85 (a), A1651 (b), A2199 (c), and those for all (d), 
early-type (e), late-type galaxies (f) of A1795 det e rmine d in §4. 21 usin g the NFW profi l e, but 



based on seve ral mass profiles: ISariderson et al.l (l2003l . thick solid). iDemarco et al.l (12003 



thick dotted! . iReiprich fc Bohringerl ( 



dot-dot-dashe d), and 



Vikhlinin et al 



20021. t hick dot-dashed). iMarkevitch et al.l Jl999l . thick 



(120061 . thick dashed). The errors of VAPs using the 
mass profile of ISanderson et al.l (120031 ) are shown by thin s olid lines. The vertical dot-da shed 
and solid lines indicate the outer significance radius, rx, in lReiprich fc Bohringerl (120021 ) and 
r2oo radius computed in this study, respectively. 
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Fig. 20.— VAPs for all galaxies of A85 (a), A1651 (b), A2199 (c), and those for all (d), 
early-type (e), late-type galaxies (f) of A1795 determined in §4.2l using the NFW profile, but 
based on several galaxy samples: galaxies with 6 < 2.0 (thick solid), galaxies with 6 < 1.8 
(thick dashed), and galaxies with 6 > (thick dotted). The errors of VAPs based on the 
galaxies with 5 < 2.0 are shown by thin solid lines. The vertical dot - dashe d and solid lines 
indicate the outer significance radius, rx, in 



Reiprich &: Bohringerl (120021 ) and r2oo radius 



computed in this study, respectively. 
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Fig. 21. — (a) VDPs and (b) projected galaxy number density profiles for the clusters in 
Sample A; (c,d) those in Sample B; and (e,f) those in Sample C. Filled circles, crosses, and 
open circles denote all, early- type, and late- type galaxies, respectively. Typical errors of 
dispersion profiles are shown by errorbars according to the subsamples. 



